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Abstract 

We determine the parton distribution functions (PDFs) in a next-to-next-to-leading order (NNLO) QCD- 
analysis of the inclusive neutral-current deep-inelastic-scattering (DIS) world data combined with the neutri- 
no-nucleon DIS di-muon data and the fixed-target Drell-Yan data. The PDF-evolution is performed in the 
Nf = 3 fixed-flavor scheme and supplementary sets of PDFs in the 4- and 5-flavor schemes are derived 
from the results in the 3-flavor scheme using matching conditions. The charm-quark DIS contribution 
is calculated in a general-mass variable-flavor-number (GMVFN) scheme interpolating between the zero- 
mass 4-flavor scheme at asymptotically large values of momentum transfer Q 2 and the 3-flavor scheme 
at the value of Q 2 = m 2 in a prescription of Buza-Matiounine-Smith-van Neerven (BMSN). The results 
in the GMVFN scheme are compared with those of the fixed-flavor scheme and other prescriptions used 
in global fits of PDFs. The strong coupling constant is measured at an accuracy of * 1.5%. We obtain 
at NNLO a s (M 2 ) = 0.1 135 ± 0.0014 in the fixed-flavor scheme and a s (M 2 ) = 0.1 129 ± 0.0014 applying 
the BMSN prescription. The implications for important standard candle and hard scattering processes at 
hadron colliders are illustrated. Predictions for cross sections of W ± - and Z-boson, the top-quark pair- and 
Higgs-boson production at the Tevatron and the LHC based on the 5-flavor PDFs of the present analysis are 
provided. 
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1 Introduction 



For many hard processes at high energies heavy flavor production forms a significant part of the 
scattering cross section. As it is well known, the scaling violations are different in the massive 
and massless cases. Therefore, in all precision measurements, a detailed treatment of the heavy 
flavor contributions is required. This applies, in particular to the extraction of the twist-2 parton 
distribution functions (PDFs) in deep-inelastic scattering (DIS). In this process 0(25%) of the 
inclusive cross section in the range of small values of x is due to the production of charm-quarks as 
measured by the HERA experiments HI and ZEUS (BO- To perform a consistent QCD-analysis 
of the DIS world data and other hard scattering data, a next-to-next-to-leading order (NNLO) 
analysis is required, which includes the 3-loop anomalous dimensions (3]| and the corresponding 
Wilson coefficients JH, in particular those for the heavy flavor contributions. The latter are known 
at leading order (LO) (5H61I and next-to-leading order (NLO) 0. In the present paper we restrict the 
analysis to the NLO heavy flavor corrections. Very recently a series of Mellin moments at NNLO 
has been calculated in Ref. (8]| for the heavy flavor Wilson coefficients of the structure function 
Fi, in the region Q 2 > 10 • mr, where is the heavy quark mass and Q 2 is the momentum transfer 
squared. Because of the large heavy flavor contribution to F2, its correct description is essential in 
precision measurements of the strong coupling constant a s and of the PDFs. 

At asymptotically large values of Q 2 , the heavy flavor contributions rise like a s (Q 2 )ln(Q 2 /m 2 ). 
Despite the suppression due to the relatively small value of a s at large scales, these terms might 
dominate and therefore their resummation is necessary [6]|. It can be easily performed through the 
renormalization group equations for mass factorization for the process independent contributions 
defining the so-called variable-flavor-number (VFN) scheme. Thereby heavy quark PDFs are in- 
troduced, as e.g. suggested in Ref. O. A VFN scheme has to be used in global fits of hadron 
collider data if the cross sections of the corresponding processes are not available in the 3-flavor 
scheme. However, since VFN schemes are only applicable at asymptotically large momentum 
transfers, one has to find a description suitable for lower virtualities Q 2 , which matches with the 
3-flavor scheme at the scale Q 2 = m 2 r , cf. Ref. IfTOl . 

At the same time, the resummed large logarithms occur in the higher order corrections. In 
the NLO corrections to the massive electro-production coefficient functions [7] the terms up to 
o: 2 (Q 2 )^ 2 (Q 2 /m 2 ) are manifest. Therefore the resummation of the remaining large logarithms is 
much less important as compared to the LO case. Furthermore, in most of the kinematic domain of 
the DIS experiments the impact of the resummation is insignificant IfTTTl . Eventually, the relevance 
of the resummation is defined by the precision of the analyzed data and has to be checked in the 
respective cases. In this paper we study the impact of the heavy flavor corrections on the PDFs 
extracted from global fits including the most recent neutral-current DIS data. We apply the results 
of the QCD-analysis to main NNLO hard scattering cross sections, as the W/Z-gauge boson, top- 
quark pair and Higgs-boson production at hadron colliders. 

The paper is organized as follows. In Section 2 we outline the theoretical formalism which 
describes the heavy quark contributions to DIS structure functions and the formulation of VFN 
schemes, cf. Refs. [8l[l2l[T3]|. A phenomenological comparison of the fixed-flavor number (FFN) 
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scheme and different VFN schemes is performed in Section 3 and 4. In Section 4 we present the 
results of an NNLO PDF-fit to the DIS world data, the fixed-target Drell-Yan- and di-muon data in 
different schemes using correlated errors to determine the PDF-parameters and a s (M^). Precision 
predictions of PDFs are very essential for all measurements at hadron colliders [fT4|. Section 5 
describes the 3-, 4-, and 5- flavor PDFs generated from the results of our fit and applications to 
hadron collider phenomenology, such as the scattering cross section of W ± - and Z-boson produc- 
tion, the top-quark pair and Higgs-boson cross sections based on the 5-flavor PDFs obtained in the 
present analysis. Section 6 contains the conclusions. 



2 Heavy Quark Contributions: Theoretical Framework 

In inclusive DIS, heavy quarks contribute to the final state if we consider extrinsic heavy flavor 
production only [JJ. In fixed-order calculations of the inclusive heavy flavor cross sections in the 
fixed flavor number scheme (FFNS) for Nf light quarks, one obtains the following representation 
for the DIS structure functions to NLO in case of single photon exchange J51|7l[8l 

F*' exact (N f ,x,Q 2 ) = 

^™\z{el\H g ,(z,Q 2 ,ml,f)-G(N 

+ ^ e 2 L g4 (z, Q 2 ,m 2 h ,f?)-G(N f , -,A + L™(z, Q 2 ,m 2 ,p 2 )-f(N f , -,p 2 )\ , 

k=l ) 

(1) 

where i = 2,L. The functions H g ^ q )j and L g ( q ),i denote the massive Wilson coefficients with the 
photon coupling to the heavy (H) or a light (L) quark line, respectively, x = Q 2 /(2p.q) is the 
Bjorken scaling variable, with q the 4-momentum transfer, p the nucleon momentum, Q 2 = -q 2 ; 
*max = Q 2 KQ 2 + ^ m j l ) is production threshold; e/ ; is the charge of the heavy quark, with h = c,b. 
We introduced a second symbol for the number of the light flavors, Ni, which counts the number 
of the light quark anti-quark final state pairs associated to the Wilson coefficients L g j. The flavor 
singlet and non-singlet distributions are given by 

Nf 

^(N f ,x,n 2 ) = Y\vk{Nf,x,ii 2 ) + q k {N f ,x,f)}, (2) 
k=l 

N f 

f(N f ,x,n 2 ) = Y j e 2 k [q k (Nf,x,i/) + q k (N f ,x,iJ 2 )] , (3) 
k=l 

where q k , q k and G are the light quark, anti-quark and gluon distributions. Here and in the following 
we identify the factorization and renormalization scales by jx = fj.f = pr. In open heavy flavor 
production one usually chooses p? = Q 2 + 4m 2 , while for the inclusive structure functions one sets 

n 2 = Q 2 . 

'Potential contributions due to intrinsic charm were limited to be less than 1 % in Ref. 03). 
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The massive Wilson coefficients in Eq. (OQ) are available in analytic form at LO [0 and in 
semi- analytic form at NLO Q q For Q 2 /m^ :» 1 they were given in analytic form to NLO in 
Refs. [fT2l[T7l[T8ll and in MM to NNLO for F L and F 2 . The NNLO contributions to F 2 are not yet 
fully available as general expressions in x or the Mellin variable N, since for one part, only a series 
of Mellin moments at fixed integer values of N has been calculated so far [8 ]. In the limit Q 2 » m 2 , 
the integration in Eq. (OQ) extends to x max = 1 and additional soft- and virtual terms contribute to 
the cross section according to the Kinoshita-Lee-Nauenberg theorem, cf. e.g. IfTTTl . 

In Ref. the effects due to heavy quark loops in external gluon lines were absorbed for the 
heavy flavor Wilson coefficients into the strong coupling constant to NLO, which is then to be 
taken in the corresponding momentum subtraction scheme in Ref. |0. The necessary changes 
for a s in the MS-scheme are discussed in Refs. j8l[l2l[T3]|. In the present paper, we will include 
the NLO contributions for Fi and F 2 with a s in the MS-scheme, cf. Ref. 10. The choice of a 
MOM scheme always forms an intermediate step, since it applies to the heavy degrees of freedom 
only. The structure functions also contain the light flavor PDFs and massless Wilson coefficients, 
the scaling violations of which are governed by only. Also, one cannot choose a scheme, 
which introduces heavy quark mass effects in the strong coupling constant below any heavy flavor 
threshold. 

In the asymptotic region Q 2 » m 2 the Wilson coefficients L g ( q ),2 and Hg( 9 ),2 for the heavy flavor 
structure function /7^> exact f gq (TJJ) can b e expressed in terms of the massive operator matrix 
elements A,-, and the massless Wilson coefficients C^2- The former are given by 
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cf. Refs. 112 [II EM 111. To NLO the massive OMEs do not depend on N f . The massless 
Wilson coefficients for the structure function F2 are given by 



C k , 2 ^Nf,z,^j = ^a n s (N f ,M 2 )C^N f ,z,^j, k = q,g, 



(7) 



cf. Refs. [31I2T1 . In case of C q p_ we decompose the Wilson coefficients into flavor non-singlet 
(NS) and pure-singlet (PS) contributions and C^|. We use the strong coupling constant in the 

notation a s {Nf,jj 2 ) = a s (Nf,fu 2 )/(47r). At the different heavy flavor thresholds /j 2 = m?, h = c,b, 
matching conditions are employed to a s {jx 2 ), cf. e.g. Ref. [|22l . 

Up to 0{a 2 ) the asymptotic expressions for the heavy flavor coefficients £g( ? ),2 and #g( ? ),2 



2 A fast implementation in Mellin space is given in Ref. |[T6l . 
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read HO 
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The symbol <g> denotes the Mellin convolution 
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(12) 



and all arguments except of Nf are omitted for brevity. Note that nearly identical graphs contribute 

to £ as y mp and the second last term of H as ^ mp . These are accounted for in different classes due to 

g,2 g,2 

the final state fermion pair, which consists of the light quarks in the first case and the heavy quark 
in the second case. Therefore we introduced Ni as a second label for the number of light flavors in 
the final state, cf. Eq. (OQ). 

The OMEs enter in the matching conditions for the PDFs in the Af /--flavor scheme with the 
ones for (Nf + 1) massless flavors 021 which are implied by the renormalization group equations. 
In particular, the NNLO heavy-quark distribution in the (Nf + l)-flavor scheme at 0(a 2 ) reads 



h w (x,fi 2 ) + h m (x,iu 2 ) = a s (N f + l,M 2 )^A[ l ^)®G (2) (Nf,fi 2 ) 



(x), 



(13) 



h {2 \x,n 2 ) + h (2) (x,n 2 ) = h w (x,fi 2 ) + h w (x,//) 



+a 2 (N f + l,iu 2 )- 



hg 
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W 2 ) 



®G (2) (N f ,t?) 



(x) + 



L (2),PS 
hq 



®2 (2) (%,^ 2 ) 



W , (14) 



where and S^ 2 ^ are the gluon and flavor singlet distributions, respectively, evolved at NNLO. 
Likewise, one obtains for the gluon, flavor non-singlet and singlet distributions in the Nf + 1 -flavor 
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scheme up to 0(a s ) 

G i2 \N f + l,x,M 2 ) = G^(N f ,x,fi 2 ) + a s (N f + \,fi 2 ) 



+a 2 s (N f + l,n 2 ) 
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These distributions obey momentum conservation 

f 

/o 



1 = 



-I 



dx x 



(16) 



and the light quark and anti-quark distributions are given by 
qf (N f + \,x,m 2 ) + q (2) (N f + l,x, M 2 ) 

®[q (2) (N f ,x,f?) + qf\Nf,x,fi 2 )] . (17) 



dxx [G(.N f ,/j 2 ,x) + X(N f ,iJ 2 ,x)] 

N f 

G(N f + 1 ,fi 2 ,x) + ^ [q k (N f + 1 , x,/u 2 ) + q k (N f +l, x, M 2 )} + h (2 \x,/u 2 ) + h (2 \x,/u 2 ) L . 



k=\ 



(18) 



Since the OMEs are process independent quantities this property is maintained by the (Nf + 1)- 
flavor PDFs. One may apply these PDFs in a hard scattering process for large enough scales 
fij, :» m^, where the power corrections are negligible. In particular, the heavy flavor structure 
function F2 is defined in the (Nf + l)-flavor scheme as the convolution of the (Nf + l)-flavor PDFs 
with the massless Wilson coefficients C q ( g ^2- This representation is the so-called zero mass VFN 
(ZMVFN) scheme expression, which is applicable only in the asymptotic region. 
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The heavy flavor part of F2 in the region Q 2 » m 2 for Nf + l flavors up to 0{a 2 s ) is given by 



./?,ZMVFN 



(N f + l,x,Q 2 ) = xel^h (2 \x,i/) + hV\x,f/) 



+a s (N f + l,fi 2 ) 



+a 2 (N f + l,ti 2 ) 



(x) 



+a s (N f + l,/f) 
Ni 



Q- 



(x) 



+xa 2 (Nf + l,fi 2 ) 



(x) 

2 1 



(X). 



(19) 



At this point we would briefly like to comment on the longitudinal structure function F^. As 
a matter of fact, the above concept of a ZMVFN scheme cannot be directly applied to the heavy 
flavor component of Fi even in the asymptotic regime of Q 2 » m?; e.g. at 0(a s ), similarly to Eq. 
(IT9T). one obtains 



^/i,asymp 



(Nf + l,x,Q z ) = a s (Nf + l,Li 2 )el 



) j Q A 



C g ,L\^)® G ( N f^) 



(X), 



(20) 



Here, the gluon density is convoluted with the LO gluon Wilson coefficient C^ L but not a splitting 
function, because unlike the case of F2 no collinear logarithm emerges. The example illustrates 
that a detailed renormalization group analysis is a necessary prerequisite to the use of heavy quark 
densities even in the asymptotic region, cf. Ref. lfT2l . 



3 Comparison of the 3- and the 4-Flavor Schemes 

At 0{a[) the universal contribution (referring to the massive OMEs only) to the heavy flavor singlet 
contribution to F2 is given by 

F^' {l \N f = 4,x,Q 2 ) = elx[h (l \x,i/) + h (l \x,iu 2 )] , 1=1,2 . (21) 

— h (1) ? ? (10) 1 — 1 7 7 

It vanishes for F 2 at /r = m h , since a hg = 0, cf. Eq. fl4J), and it is negative for jx < m h . 
However, the 1st order heavy quark contribution to the structure function F2 is positive, since 
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Figure 1: Matching of F^ BMSN {N f = 4,x,Q 2 ) (solid lines) with F c 2 exact (N f = 3,x,Q 2 ) (dash-dotted lines) 
at small Q 2 in 0(a 2 ). The vertical line denotes the position of the charm-quark mass m c - 1.43 GeV. 

the /j 2 dependence is canceled by a corresponding logarithm oc \n(Q 2 //j. 2 ) in the massless Wilson 
coefficient C^l in Eq.(fT9l). Despite that in the 3-fiavor scheme the heavy quark contribution to 
the DIS structure functions also falls at small Q 2 , it is present down to the photo-production limit. 
At 0(a 2 ) the agreement between the two schemes at low Q 2 is even worse since the term a hg is 

— h (2) 9 

negative which implies F. (Nf = 4,x, Q ) < because the gluon contribution dominates over the 
pure-singlet part numerically. The impact of the large-log resummation is negligible at small scales 
Q 2 and any reasonable scheme must reproduce the 3-fiavor scheme. Therefore, at low values of 
Q 2 the ZMVFN scheme is not applicable. It has to be modified according to practical purposes. 
VFN schemes with such modifications are called general-mass variable-flavor-number (GMVFN) 
schemes, in contrast to the ZMVFN scheme. A particular form of the GMVFN scheme cannot 
be derived from first principles in an unique way, but is subject to the corresponding prescription. 
Consistent schemes have to obey renormalization group equations to not violate the running of the 
coupling constant and masses, and to obey correct scale evolution. As a general requirement any 
such prescription should provide a continuous transition from the 3-fiavor scheme at low values of 
/j 2 to the 4-fiavor scheme at large scales. 

An early formulation of a GMVFN scheme by Aivazis-Collins-Olness-Tung (ACOT) [9] does 
not allow a smooth matching with the 3-flavor scheme at small scales Q 2 . In the ACOT scheme, the 
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Figure 2: Comparison of F c 2 in different schemes to HI- and ZEUS -data. Solid lines: GMVFN scheme 
in the BMSN prescription, dash-dotted lines: 3-flavor scheme, dashed lines: 4-flavor scheme. The vertical 
dotted line denotes the position of the charm-quark mass m c = 1.43 GeV. 



slope in Q 2 turns out to be too large. Later the so-called Thorne-Roberts (TR) scheme overcoming 
this shortcoming was suggested [23]. However, beyond NLO this scheme is very involved and its 
numerical implementation is problematic It2"4l . Recently, the early ACOT prescription has been 
modified in order to improve the behavior at low values of Q 2 [25]. This modified description, 
the so-called ACOTCf) scheme, is used in particular at NNLO in Ref. [J26]| . Another GMVFN 
prescription, which was suggested earlier by Buza-Matiounine-Smith-van Neerven (BMSN) |[T2l 
for F 2 , is defined by 

F h 2 BMSN (N f + l,x,Q 2 ) = F h *™\N f ,x, Q 2 ) + Ff MVFN (V/ + l,x, Q 2 )- F h 2 ^(N f ,x, Q 2 ) , (22) 
with Nf = 3 for h = c. 

Note that the difference of the last two terms in Eq. (1221 depends on Nf through the strong 
coupling constant only, which is a specific feature up to NLO. For the choice of /u 2 = Q 2 the 
asymptotic terms cancel at Q 2 = m 2 h in Eq. (1221 . In this limit F 2 ' BMSN (Nf = 4) reproduces the result 
in the 3-flavor scheme. Moreover, F 2 ' B (Nf = 4) matches with the 3-flavor scheme smoothly 
as shown in Figure HJ Minor kinks between F^ BMSN (N f = 4) and F 2 ' sx<iCt (Nf = 3) stem from the 
matching of a s (Nf,/u 2 ) at jx 2 = m 2 . It appears since the matching condition for a s (Nf,jx 2 ) does 
provide a continuous but not a smooth transition at the flavor thresholds. The numerical impact 
of this kink is marginal in the analysis of the current data. At large Q 2 the asymptotic expression 
F 2 ' asymp (N f = 3) cancels the term F 2 ' e ™ ct (N f = 3) in Eq. d22j and Fj™ 8 ^/ = 4) reproduces 
the result in the ZMVFN scheme. The cancellation is not perfect due to the difference in the 
upper limit of integration in Eq. (QQ) and the expression for the ZMVFN scheme, which affects 
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only the non-singlet Compton-type contribution given by the coefficient functions L^. For the 
3-flavor expression of Eq. this term rises as ln 3 (<2 2 /m?) at large Q 2 . In the asymptotic limit 
of Ref. ifm the corresponding singular contribution is washed out in ^ s y m P ,NS As a result, there 

remains a contribution ~ ln 3 (g 2 /m 2 ) in the difference of F% exact (N f = 3) and F^ symp (N f = 3). 
This mismatch is caused by a well-known soft and virtual term, which occurs in the inclusive 
analysis for large arguments of the Wilson coefficient and is easily corrected, cf. Ref. [117112711 . On 
the other hand, the non-singlet contribution to heavy quark electro-production is numerically very 
small and the term ~ ln 3 (<2 2 /m 2 ) is apparent only at very large values of Q 2 and relatively large x. 
The accuracy of realistic data at this kinematics is rather poor and even for the definition of Eq. CD 
the impact of the mismatch between F^ e ™ ci (Nf = 3) and F 2 ,asymp (iV/ = 3) turns out to be marginal 
in the data analysis. 



x=0.002 




_ _ 1 10 10 

Q (GeV 2 ) Q (GeV 2 ) 



Figure 3: The c-quark distributions calculated using the fixed-order relation Eqs. (1 141 13b (solid lines) 
compared to the result in the 4-flavor scheme evolving from m 2 c and using Eqs. (1 14113b as boundary condition 
(dash-dotted lines) at 0(a s ) (left panel) and at 0(a 2 ) (right panel). 



A representative set of the ZEUS- and Hl-data iH on F c 2 is compared to F^ (Nf = 4), 
F% exact (N f = 3), and F% ZMWFN (N f = 4) in Figured The 3-flavor PDFs used in this comparison are 
evolved starting from m c = 1.43 GeV with the input given by the MRST2001 PDFs of Ref. BH. 
Because of the kinematic constraints, at x ~ 0.0001 only the values of Q 2 < 10 GeV 2 are available 
in the data. At such low scales the calculation in the 3-flavor and BMSN scheme yield practically 
the same results. At x ~ 0.01 the typical values of Q 2 are much bigger. In this kinematic region, 
the BMSN scheme yields a larger contribution than obtained in the 3-flavor scheme. However, 
the uncertainties in the data are still quite large due to limited statistics. The comparison of the 
calculations with the data is rather insensitive to the choice of scheme. The non-singlet term in 
Eq. (OQ) is not taken into account in the comparisons shown in Figure[2l Its impact is most significant 
at large values of x and Q 2 , but even in this case, it is much smaller than the data uncertainty. For 
intermediate values of x ~ 0.001, a combination of these two cases is observed: at large Q 2 the 
uncertainties in the data do not allow to distinguish between both schemes, while at small Q 2 
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Figure 4: Comparison of predictions in different schemes to ZEUS- and Hl-data on F^{x, Q 2 ). The 
notations are the same as in Figure |2] The vertical dotted line marks the position of my - 4.3 GeV. 



the numerical difference between the 3-flavor and the BMSN scheme calculations is very small. 
Summarizing, for the analysis of realistic data on F c 2 , the BMSN scheme is very similar to the 
3-flavor scheme. This is not a particular feature of the BMSN prescription, since the difference 
between 3- and 4-flavor schemes at large Q 2 is also smaller than the uncertainties in the available 
data and, once the smooth matching is provided, a GMVFN scheme must be close to the 3-flavor 
one at small Q 2 . This conclusion is in agreement with the results of Ref. 11301 . It derives from the 
fact that once the 0(a 2 ) corrections are taken into account the need of a large-log resummation 
is thus greatly reduced, which is well known for a long time, cf. Ref. [11]. In Figure [3] the c- 
quark distribution defined in Eq. (1131141) is compared to the one evolved in the 4-flavor scheme 
starting from the scale of m c using Eqs. (1131141) as boundary conditions. The former is derived 
from fixed-order perturbation theory, while for the latter resummation is performed through the 
evolution equations. At 0(a s ) the difference between these two approaches is significant indeed, 
however, at 0(a 2 ) it is much smaller, and quite unimportant for realistic kinematics. 

As evident from Figure [2] the scheme choice cannot resolve observed discrepancies between 
data and the theoretical predictions to NLO. Given the mass of the charm quark and the PDFs 
determined in inclusive analyses, higher order QCD corrections are needed. In particular, at small 
x and Q 2 , the partial 0{al) corrections to the massive Wilson coefficient H g 2 obtained through 
threshold resummation [[3"T]| give a significant contribution to F c 2 and greatly improve the agreement 
to the data 11321 . In this kinematic region, the integral of Eq. <Q]) is mostly sensitive to the threshold 
of heavy quark production and the approximate form of H g 2 derived in Ref. 1(321 is sufficient. At 
large values of Q 2 , the threshold approximation is inapplicable, and a complete NNLO calculation 
is required, cf. Ref. (S]|. For &-quark production the resummation effects are less important since 
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Figure 5: The errors in the inclusive structure function Fiix, Q 1 ) measured by the HI collaboration ll33l 
in comparison with the impact of the heavy quark scheme variation on the QCD calculations for F2&, Q 2 ). 
Solid line: c-quark contribution, dash-dotted lines: fr-quark contribution. The vertical dots mark the posi- 
tions of m c = 1.43 GeV and nib = 4.3 GeV, respectively. 

the asymptotic region is scaled to bigger values of Q 2 and the data are less precise due to the 
smaller scattering cross section. This is illustrated by a comparison of the ZEUS data on with 
calculations of the 4-flavor ZMVFN scheme, the 3-flavor scheme, and the BMSN prescription for 
the GMVFN scheme given in Figure [4] 

Also the inclusive structure function F2 is sensitive to the choice of the heavy quark scheme, 
due to the significant charm contribution in the small x region. In fact, for F2 the sensitivity is 
much larger than for the heavy quark contributions F c 2 and F^ alone due to the far higher accuracy 
of the data. In Figure [5] we compare the errors in F2 measured by the HI collaboration [33] with 
the difference between p^' exact _ p > For £-quark production the scheme variation effect, 

which is calculated in the same way as in Figure HI is negligible as compared to the accuracy of 
the data in the entire phase space. For the c-quark contributions, maximal sensitivity to the scheme 
choice appears at largest values of Q 2 at x ~ 0.001, similarly to the case of the data for F c 2 given 
in Figure [21 The effect is localized in phase space and appears to be at the margin of the statistical 
resolution. Therefore the impact of the scheme variation on the data analysis turns out to be rather 
mild. To check it in a more quantitative way we compare the QCD-analysis of the inclusive DIS 
data performed in the 3-flavor scheme with the one in the BMSN prescription of the GMVFN 
scheme. Details and results of these analyses are described in the following Section. 
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4 Impact of the Scheme Choice on the PDFs 



We determine the PDFs from the inclusive DIS world data obtained at the HERA collider and 
in the fixed-target experiments ll33l[34ll . These data are supplemented by the fixed-target Drell- 
Yan data ll35l and the di-muon data from (anti)neutrino-nucleon DIS 11361 . which allow the flavor 
separation of the sea quark distributions. Details of the data selection, the corrections applied to 
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71 


72 


A 




0.662 + 0.034 


3.574 + 0.078 


-0.590 + 0.027 


-0.71+0.17 




d v 


1.06 + 0.12 


6.42 + 0.41 


4.4 + 1.0 


-7.0+1.3 




u s 


-0.216 + 0.011 


6.83+0.24 


0.64 + 0.29 




0.1408 + 0.0079 


A 


0.7 


11.7 + 1.9 


-3.5 + 2.1 




0.256 + 0.082 


s 


-0.253 + 0.058 


7.61+0.65 






0.080 + 0.016 


G 


-0.214 + 0.013 


7.95 + 0.15 


0.65 + 0.92 







Table 1: The parameters of the PDFs and their lcr errors in the 3-flavor scheme. 



the data, and statistical procedures used in the analysis can be found in Refs. H371I38I . The analysis 
is performed taking into account the NNLO corrections for the light flavor contributions. 

For the neutral-current c-quark contributions to the structure function F% two variants are com- 
pared, both up to the level of the 0{a 2 ) corrections, d In one case we employ i 7 ^' exact (./V/ = 3) of 
Eq. (OQ), calculated for three light quark flavors choosing the factorization scale of /j 2 = Q 2 + 4m 2 .. 
This is compared to the BMSN prescription of the GMVFN scheme F c ^ BMSN (N f = 4) given by 
Eq. (1221) with the factorization scale /j 2 = Q 2 . However, in the kinematical region of the data this 
variation of scale yields no difference in the fit results. Our fit is based on the reduced cross sec- 
tions rather than the DIS structure functions. Therefore we also have to consider the longitudinal 
structure function Fl- Since the data are much less sensitive to Ft than to F2 the scheme choice 
is unimportant for the former and in both variants of the fit it is calculated in the Nf = 3 FFN 
scheme, Eq. (OQ)- Likewise this is the case for the £-quark contribution, where the scheme choice 
is also unimportant, as one can see in the comparisons of Section [3j The Nf = 3 FFN scheme 
is used both for F h 2 and F h L . The charged-current c-quark contribution to the structure functions, 
which are related to the di-muon (anti)neutrino-nucleon DIS data used in the fit, are calculated in 
the Nf = 3 FFN scheme at NLO 11391 . For the Drell-Yan cross sections, we include the NNLO 
QCD corrections Il40l . In this case the 5-flavor PDFs defined in Eqs. (fT4T - fT71) are used in order to 
take into account the c- and Z?-quark contributions. Note, however, that at the typical fixed-target 
energies the impact of heavy quarks is marginal and the 3-flavor scheme provides a sufficiently 
good description. 

The proton PDFs are parameterized at the scale Qq = 9 GeV 2 in the 3-flavor scheme. At 
the starting scale the following functions are used for the valence quark, gluon, and sea-quark 

3 The effects of the O(a^) corrections calculated recently in Ref. [8 j will be studied in a forthcoming paper. 
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distributions : 



xqv(x,Q 2 ) = qi ' v - x a "(l-x) b "x p "' vix \ P g ,v = 7h q x + y 2 , q x 2 , q = u,d, 



N v 



(23) 



xG(x,Q 2 ) = A G x ac (l-x) bc x PG{x \ P G = 7l,GX, (24) 

xu s (x,Q 2 ) = xu s {x,Q 2 )=A u x a "%l-x) b " s x p ^ {x \ P u , s =7l,u S x, (25) 

xA(x,Q 2 ) = xd s (x,Q 2 )-xu s (x,Q 2 )=A A x a \l-x) b *x p ^ x \ P A = y UA x. (26) 

The strange quark distribution is taken in the charge-symmetric form 

xs(x, Q 2 ) = xs(x, Gq) = A^'Cl - x)\ (27) 



in agreement with the results of Ref. 11381 . The polynomials P(x) used in Eqs. (12314261) provide 
sufficient flexibility of the PDF-parameterization with respect to the analyzed data and no addi- 
tional terms are required to improve the fit quality. The PDF-parameters determined from the fit 
performed in the 3-flavor scheme are given in Tabled! Because of the lack of the neutron-target 
data in the region of small values of x, the low-* exponent a A cannot be defined from the fit and 
we fix it to 0.7 to choose an ansatz, in agreement with the values obtained for the low-x exponents 
of the valence quark distributions and phenomenological estimates, cf. e.g. |]4Tfl. However, once 
we have fixed a A the uncertainty in the sea quark distributions at small x is underestimated. We 
therefore choose an uncertainty 6a A = 0.3, and, in order to account for its impact on the other 
PDF-parameters, we calculate the errors in the latter with the value of a A released, but with an 
additional pseudo-measurement of a A = 0.7 + 0.3 added to the data set. In our fit the heavy quark 
masses are fixed at m c = 1.5 GeV and m\, = 4.5 GeV and the same approach is employed to take 
into account possible variations of m c and m& in the ranges of +0.1 GeV and +0.5 GeV, respec- 
tively. Note that the normalization parameters for the valence quarks and gluons are defined from 
other PDF-parameters applying both fermion number- and momentum conservation. In the global 
fit we obtain 

X 2 3038 

-^— = = 1.1 (28) 

NDP 2716 

for the parameter values listed in Table [T] 

In the fit 25 parameters are determined. The covariance matrix elements for these parameters 
are given in Table 2. The parameter errors quoted are due to the propagation of the statistical and 
systematic errors in the data. The error correlations are taken into account if available, which is 
the case for most of the data sets considered. 

The gluon and flavor singlet distributions obtained in case of the BMSN prescription are com- 
pared to those referring to the 3-flavor scheme in Figure[6] The difference between the two variants 
is quite small and situated well within the PDF-uncertainties. For the non-singlet PDFs it is even 
smaller, since the heavy quark contribution is negligible at x > 0.1, cf. Ref. H21 . For the BMSN 
variant of the fit a value of ^ 2 /NDP = 3036/2716 is obtained, very close to the one for the fit in 
the 3-flavor scheme. This is in line with the comparisons given in Section [31 which show that in 
the case of a smooth matching of the 3-flavor and VFN scheme at small values of Q 2 , there is little 
room for a difference between them in the region of the present experiments. 
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(1=3 GeV 




Figure 6: The lcr error band for the gluon (left panel) and sea (right panel) distributions obtained in two 
valiants of the fit. Solid lines: 3-flavor scheme, dashed lines : GMVFN scheme in the BMSN prescription. 



The difference between the fits performed using the TR-prescription for the GMVFN scheme 
and in the 3-flavor scheme is also not dramatic, as one can see in the ZEUS NLO PDF fit, Ref. Il43l . 
However, it is somewhat larger than in the case of the BMSN prescription. In particular, this 
happens since the TR-prescription does not provide a smooth matching with the 3-flavor scheme at 
low values of Q 2 and at small values of x. By construction, the TR-prescription provides a smooth 
transition for the gluon-initiated contribution only. However, at small values of Q 2 the gluon 
distribution has a valence-like form and falls at small x. As a result, the quark-singlet contribution 
to the slope dF c 2 l d\n(Q 2 ) is non-negligible at small values of x, which leads to a kink at the 
matching point Q 2 = m 2 h in F\ using the TR-prescription, see Figure [7] This is an artefact of the 
description leading to an overestimation of the heavy quark contribution and, correspondingly, an 
underestimation of the fitted light quark PDFs at small values of x. In the ACOT(x)-prescription the 
smoothness is not required by definition and the kink in is even bigger than for the case of TR- 
prescription. In the most recent version of the ACOT-prescription this problem is addressed B4l 
and should yield a result closer to the 3-flavor scheme than the ACOT(^)-prescription. 

Summarizing the comparisons of Sections[3]and|4]we conclude that, once the 0{a 2 ) corrections 
to heavy quark electro-production are taken into account and a smooth matching with the 3-flavor 
scheme at small Q 2 is provided, the GMVFN scheme should agree to the 3-flavor scheme for the 
kinematics explored by experiments so far. Furthermore, it is expected that the NNLO corrections 
to the heavy quark structure functions of Ref. [8] lead to an even better agreement. 
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x=0.0001 



0.12 



3-flavour 



Thorne-Roberts 




Figure 7: Matching of F% TR (x, Q 2 ) (dash-dotted line) and F' 2 

3,x,Q 2 ) (solid line) at small Q 2 at 0(a s ). The MRST2001 PDFs of Ref. (H are used. The vertical line 
denotes the position of the charm-quark mass m c = 1.43 GeV 



r c,ACOTCtr) 



Q 2 (GeV 2 ) 

(x, Q 2 ) (dashed line) with F c 7 ' exact (N f = 



For the strong coupling constant at NNLO in QCD the values 

af s (N f = 5,M z ) = 0.1 135 + 0.0014 (exp) FFN scheme, N f = 3 (29) 
af s (N f = 5,Mz) = 0.1 129 + 0.0014 (exp) BMSNscheme (30) 

are obtained. The small difference between these two values lies well within the experimental 
uncertainty. In Table [3] we compare these values to other recent NNLO determinations of the 
strong coupling constant. Our results agree very well with those of Refs. [|42ll45ll . Note that the 
data sets used in the non-singlet fit of Ref. Il42l are rather different from those used in the present 
analysis. The value of or 5 (M|) given in Ref. [|46l is by 2.1a larger. As is well-known from the 
non-singlet data analysis 11421 . a somewhat higher value of a s (M^) is obtained at N 3 LO, cf. also 
Ref. [21] for an estimate. The difference of these determinations at NNLO and N 3 LO is half of 
the experimental error found in the present analysis. Eqs. d29l and (l30l) determine a s (M^) at an 
accuracy of ~ 1.5%. 
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Figure 8: Left panel: The 5-flavor PDFs at the scale of /j 2 = 10 4 GeV 2 and the 3-flavor gluon distributoin 
given for comparison. Right panel: The ratio of the 5-flavor to 3-flavor gluon (solid line) and singlet (dash- 
dotted line) distributions at the same scale. 



5 Applications to Collider Phenomenology 

In this Section, we investigate the implications of the PDFs obtained in the present NNLO analysis 
for collider phenomenology. To that end we focus on important (semi)-inclusive scattering cross 
sections at hadron colliders, such as the Drell-Yan process for W ± - and Z-boson production, the 
pair-production of top-quarks and (Standard Model) Higgs-boson production. The corresponding 
cross sections in, say, proton-proton scattering can be written as 

(Tpp-+x(s) = 2^ I dxidx2fi(xi,v 2 )fj(x 2 ,iJ 2 )d-ij^x(xuX2,s,a s (ju 2 ),iu 2 ) , (31) 
ij J 

where X is the final state under consideration and s the c.m.s. energy. The PDFs are collectively 
denoted by fi, fj and the sum runs over all partons. At hadron colliders the convolution of fi and 
fj parameterizes the so-called parton luminosity L ;J . In the following, we will employ our PDFs 
and present numbers for pp-collisions at Tevatron with ^ = 1.96 TeV and for pp-collisions at 
LHC at energies ^ = 7,10 and 14 TeV. To that end, we have to rely on the perturbative QCD 
evolution of the light and heavy PDFs to Tevatron and LHC scales, which puts us also in the 
position to compare to other global PDF analyses. In the comparison we will consider the impact 
of the error of the PDFs at the level icrp (the index P denoting PDFs), which results from the 
experimental errors in the PDF analysis, including full error correlation, see Section|4j We will not 
consider theory errors implied by varying the factorization and renormalization scales. At the level 
of NNLO these amount typically to a few per cent only and, moreover, are largely independent of 
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a,(M|) 




ABKM 


0.1135 + 0.0014 


heavy quarks: FFN Nf = 3 


ABKM 


0.1129 + 0.0014 


heavy quarks: BMSN approach 


BBG Ea 


„ , „ „ , +0.0019 
0.1134 

-0.0021 


valence analysis, NNLO 


AMP06 1571 


0.1128 + 0.0015 




JRII451 


0.1124 + 0.0020 


dynamical approach 


MSTW||46ll 


0.1171 ±0.0014 




BBG BH 


0.1141 +0 - 002 ° 
-0.0022 


valence analysis, N 3 LO 



Table 3 : Comparison of different measurements of a V (M|) at NNLO and higher order. 



the PDFs. Also the anticipated statistical and systematic errors in the measurements at Tevatron 
and the corresponding resolutions, which can be achieved at the LHC, are not considered. 

5.1 Evolution of Light and Heavy PDFs 

The typical energy scales for hard scattering processes at high-energy hadron colliders are often 
much larger than the c-quark mass, and even than the &-quark mass. In this case the (4- )5-flavor 
scheme is the relevant choice, if power corrections and non-factorizing contributions can be safely 
neglected. Moreover, very often this is the only approach feasible, since the cross sections of the 
partonic sub-processes are only available in the approximation of massless initial-state partons. 
The 3-flavor PDFs obtained from the fit in Section|4]can be used to generate the 4-flavor distribu- 
tions using the matching conditions in Eqs. (fl4l4T71) . As we show in Figure [31 at 0(a 2 ) and low 
scales, the PDFs computed in this way are very similar to the evolved ones, provided the matched 
PDFs are taken as boundary conditions in the evolution. At large scales the difference between 
these two cases is non-negligible, contrary to the case of heavy quark DIS electro-production. 
The large-log resummation effects can be important in some range of the phase space at hadron 
colliders. In view of these aspects we obtain the 4-flavor PDFs from the NNLO evolution with 
the boundary scale m 2 c and the boundary conditions given in Eqs. (fT4l4T71) . The 5-fiavor PDFs 
are obtained from the evolved 4-flavor distributions using analogous boundary conditions at the 
scale m?. Since at m, the mass effects of the charm-quark are not negligible due to mym c ~ 10, 
this approach is some approximation, the validity of which has to be tested for the corresponding 
processes. The problem of the heavy quark mass scale separation cannot be resolved within the 
concept of a VFN scheme and implies an unavoidable theoretical uncertainty related to the use of 
VFN PDFs. In general, the heavy quark PDFs rise with the scale p 2 , while the (Nf + l)-light 
PDFs decrease correspondingly with respect to A/y-light PDFs. At the scale of 10 4 GeV 2 the 5- 
fiavor gluons lose some 7% of momentum as compared to the 3-flavor ones. This momentum is 
transferred to the c- and b-quark distributions, see Figure [8] The difference between the 3-flavor 
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singlet distribution ~E(Nf = 3) and the 5-flavor distribution 

3 

£/ K = 5 ) = Z K = 5 ) + * fa = 5 )1 (32) 

jfe=i 

is smaller than that for the gluons since in the quark-case the corresponding OMEs appear only at 
0(a 2 ). At small values of x the 5-flavor light quark and gluon distributions receive an additional 
enhancement as compared to the 3-flavor distributions due to evolution, see Figure [9] This differ- 
ence can be considered as an estimate of the theoretical uncertainty in the 5-flavor PDFs due to the 
higher order corrections. For the c- and b-quark distributions at x ~ 0.1 the effect of the evolution 
is much larger. However, due to the smallness of the heavy quark PDFs in this region its absolute 
magnitude is insignificant for most practical purposes. 



5.2 Comparison with Other NNLO Analyses 

In Figures [10] and Qj] we compare the NNLO PDFs obtained in the present analysis to the PDFs 
by Martin-Stirling-Thorne-Watt of 2008 (MSTW 2008), |47]. At the scales of /u 2 = 100 GeV 2 and 
/j 2 = 10 4 GeV, we compare the 5-flavor PDFs and at the scale of /j 2 = 4 GeV 2 , the 4-flavor PDFs 
since for the MSTW2008 set the number of flavors is four at m c < fi < m\, and 5 at /i > m^B At small 

4 If the scale is not much larger than m^. the choice of 3-flavor PDFs is most relevant, cf. Sections [213 1 
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Figure 10: The lcr error bands (shaded area) for our NNLO 4-flavor (left panel) and 5-flavor (central and 
right panels) u-, cf-quark, and gluon distributions in comparison to the corresponding MSTW2008 NNLO 
distributions B71 (dashed lines). 



values of x our gluon distribution is larger than that of MSTW2008. This difference is particularly 
essential at smaller scales where the NNLO MSTW2008 gluon distribution becomes negative at 
x ~ 5 ■ 10~ 5 . This is not the case in our analysis. Also our sea-quark distributions are larger than 
those of MSTW2008 in the small-* region. As we discuss in Section|4l this might be partly related 
to the heavy quark contribution in the GMVFN scheme employed in the MSTW2008 fit. The shape 
of the gluon distribution at small x is sensitive to the recent measurements of Fl at small Q 2 by 
the HI and ZEUS collaborations ||48l . These measurements are in agreement with our shape for 
the gluon and do slightly disfavor the MSTW2008 predictions. At large x the MSTW2008 gluon 
distribution is somewhat larger than ours due to the impact of the Tevatron jet data included in the 
MSTW2008 analysis. 

In Figure [12] we compare our 3-fiavor PDFs to the results obtained by the Dortmund group 
[Jimenes-Delgado and Reya (JR)] in Ref. 03. for j" 2 = 4, 100 and 10000 GeV 2 . At /u 2 = 4 GeV 2 
the gluon PDF ||45l is somewhat smaller for x < 5 • 10~ 5 than the gluon distribution determined in 
the present fit. This is a region in which the fit is not constrained by data. A very small difference 
is also observed for the u- and J-quark distributions in the region x ~ 0.1. Otherwise, one notices 
very good agreement of both distributions. 

In Table |4] we summarize different values of the 2nd moment of the valence quark densities. 

5 We thank P. Jimenez-Delgado and W.J. Stirling for providing us with the moments of the JR and MSTW08 
distributions. 

6 Here and in the following we restrict the comparison to the results obtained in NNLO analyses. Currently available 
NLO analyses (see in Ref. l42ll and Refs. 1133114914511 ) contain relatively large theory uncertainties of ±0.0050 for 
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Figure 11: The Icr error bands (shaded area) for our NNLO 4-flavor (left panel) and 5-flavor (central 
and right panels) s-, c-, and Z?-quark distributions in comparison to the corresponding MSTW2008 NNLO 
distributions ll47l (dashed lines). 



They are closely related to the moments which are currently measured in lattice simulations [[521 . 
The values of all analyses are very similar, with some differences still visible. A quantity of central 
importance is 

{xV(Q 2 )) = ^ dxx{[u(x,Q 2 ) + u(x,Q 2 )]-[d(x,Q 2 ) + d(x,Q 2 )]} . (33) 

In the present analysis we obtain 

(xV(Ql)) = 0.1646 + 0.0027 (this analysis) (34) 
(xV(Q 2 )) = 0.1610 + 0.0043 N 3 LO , (35) 

f° r Qq = 4 GeV 2 , where we combine in Eq. (T331) the value of the difference x(u v - d v ) obtained in 
Ref. 021 with the value for 

(x[d -u]) = 0.0072 + 0.0007 (36) 

found in the present analysis. In the above combination the correlation to the heavier flavor distri- 
butions is negligible. 

The PDF-uncertainties given in Figures [T0lfT2l are defined by the uncertainties in the analyzed 
data and the uncertainties due to m c , mb, and the low-x non-singlet exponent as discussed in 
Section HI For the c- and ^-quark distributions, the essential uncertainties are due to m c and m^, 

a s (M|), much larger than the experimental accuracy presently reached. 
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Figure 12: The lo~ error bands (shaded area) for our NNLO 3-flavor u-, d-, s-quark and gluon distributions 
in comparison to the corresponding JR NNLO distributions P31 (dashed lines). 

respectively. At small x, however they are determined much more precisely than the strange sea 
distribution, which is widely unconstrained at x < 0.01 by the present data. We now turn to some 
important inclusive processes at hadron colliders, for which we illustrate the impact of NNLO 
PDFs derived in the present analysis. 

5.3 W/Z-boson Production 

The inclusive production cross sections of single W ± - and Z-bosons are considered so-called stan- 
dard candle processes at hadron colliders. The cross sections and distributions for these processes 
are calculated up to the NNLO H40U53J - [55l (see also Ref. [56] for the expressions in Mellin-space) 
which allows to reduce the theoretical uncertainty due to factorization and renormalization scale 
variation down to a few per cent. With this theoretical accuracy provided the measurement of 
the W^/Z-boson production rates can be used to monitor the luminosity of the collider. More- 
over, a combination of the data on W^/Z-production with the non-resonant Drell-Yan data allows 
to separate the quark distributions of different flavors with a very good accuracy, cf. 11571 . The 
quark-anti-quark luminosities contributing to W + -production in pp collisions are given by 

Lqq = T[q{^Re Y Mw)q(^e~ Y Mw) + q(^Re Y ,Mw)q(^Re- Y ,Mw)}, (37) 

where r = M^/s, s denotes the c.m.s. collision energy squared, and Y is the W + c.m.s rapidity. 
In Figure [13] we compare the luminosities of Eq. (1371) weighted by the corresponding Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements . for different channels at the energy of the LHC 
with our NNLO 5-fiavor PDFs used as input and 
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(xdv(x)) 


(xlu„ — d^A(x)) 

L V "VJV // 


ABMK 


0.2981+0.0025 


0.1191 + 0.0023 


0.1790 + 0.0023 


BBGlH 


0.2986 ± 0.0029 


0.1239 + 0.0026 


0.1747 + 0.0039 


jrhs 


0.2900 ±0.0030 


0.1250 + 0.0050 


0.1640 + 0.0060 


mstw mi 


+0.0051 

0.2816 

-0.0042 


+0.0027 

0.1171 

-0.0028 


+0.0046 

0.1645 

-0.0034 


AMP06 El 


0.2947 + 0.0030 


0.1129 + 0.0031 


0.1820 + 0.0056 


BBG [N 3 LO] EH 


0.3006 + 0.0031 


0.1252 + 0.0027 


0.1754 + 0.0041 



Table 4: Comparison of the 2nd moment of the valence quark distributions at NNLO and N 3 LO obtained 
in different analyses at Q 2 = 4 GeV 2 . 5 



V* [TeV] 


this paper 


MSTW E71 


cr(W + +W-) 


aiZ) 


o-(W++W~) 


a(Z) 


1.96 (pp) 


26.2 + 0.3 


7.73 + 0.08 


25.4 + 0.4 


7.45 + 0.13 


1 (pp) 


98.8 + 1.5 


28.6 + 0.5 






10 (pp) 


145.6 + 2.4 


42.7 + 0.7 


142.1 + 2.4 


42.5 + 0.7 


14 (pp) 


207.4 + 3.7 


61.4+ 1.1 


201.1 + 3.3 


61.0+1.0 



Table 5: The total W ± and Z-cross sections [nb] at the Tevatron and LHC at the scale p = Mw/z (see 
Eq. (1381 ) for the other parameters) with the PDFs and its estimated uncertainties from the present analysis 
and in comparison to results of Ref. B71 . 

= V% = 0.9474, V% = 0.0509, Mw = 80.398 GeV . (38) 

In the forward region of rapidity the main contribution comes from the u-d annihilation. In the 
central region the c-quark contribution is also essential. Therefore, the single W ± cross section 
measurement can be used to check the magnitude of the c-quark distribution. For the case of 
antiproton-proton collisions, the quark- anti-quark luminosities are similar to Eq. (1371) ; however, 
at Tevatron the valence M-quark contribution is dominating in the whole range of rapidity. The 
cross sections for W^/Z-production at the scale /i = Mw/z for the parameters in Eq. (|38l ). and 
M z = 91.188 GeV, and including the NNLO corrections of Ref s. 15511551 are given in Table The 
quoted uncertainties are propagated from the uncertainties in the parameters of our PDFs, a s , m c 
and nib, cf. Section |4j They amount to ~ 1% at the Tevatron and ~ 2% at the LHC. Comparing the 
present analysis to Ref. [47] the results for Tevatron are at variance by 2crp, while the same cross 
sections are obtained for Z-boson production at LHC energies. 
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Figure 13: The lo~ band for the quark-anti-quark luminosities contributing to the W + production in the 
proton-proton collisions at the c.m.s. energy of V? = 14 TeV (left panel) and antiproton-proton collisions at 
the c.m.s. energy of -^s - 1-96 TeV (right panel). 



5.4 Top-quark Pair-Production 

The scattering cross section for hadroproduction of heavy quarks of mass m/ 7 is known exactly in 
QCD including radiative corrections at NLO ll58l - [6Tll . At NNLO approximate results based on the 

complete logarithmic dependence on the heavy quark velocity fi = J I -4m^/s near threshold s - 

4m? (s being the partonic c.m.s. energy) together with the exact dependence on the scale jx provide 
currently the best estimates [[621I631 . At Tevatron, the cross section is most sensitive to the qq- 



V* (TeV) 


this paper 


MSTW2008 


1.96 (pp) 


6.91+0.17 


7.04 


v (pp) 


131.3 + 7.5 


160.5 


10 (pp) 


343 + 15 


403 


14 (pp) 


780 + 28 


887 



Table 6: The total ??-production cross sections [pb] at the Tevatron and LHC for a pole mass of m t = 
173 GeV at the scale /i = m t . The results for the PDFs and its estimated uncertainties from the present 
analysis are compared to the central values obtained using the PDFs of Ref. BTl . 



annihilation channel, with the luminosities Ljj ordered in magnitude according to L qq > L qg > L gg . 
At the LHC, on the other hand, the cross section receives the dominant contribution from the gg- 
channel, in particular, from the gluon PDF in the region x « 2.5 • 10~ 2 . This makes the cross section 
for top-quark pair-production an interesting observable to investigate the gluon luminosity. Also 
the correlations of rates for /f-pairs with other cross sections can be studied quantitatively Il64l . 
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Our cross sections for ^-production are summarized in Table [6] for a pole mass of m t = 173 
GeV. We estimate the relative accuracy due to the PDF-fit for Tevatron by ~ 3%, and for the LHC 
by ~ 3.5 -4.5%. With comparison of the cross sections obtained with the PDFs of Ref. H71 we 
find agreement within lap for Tevatron. For LHC energies, the results for the MSTW08 set are 
larger by 4o> due to a bigger value of a s (M 2 ) and the larger value of the gluon PDF in the partonic 
threshold region s - 4m 2 . Note that the variation of the factorization and renormalization scale is 
not considered here. It contributes separately to theoretical uncertainty (at NNLO ~ 3 - 4% at 
Tevatron and LHC, see H62U63B for details). 

5.5 Higgs Boson Production 

Higgs-boson production is the most prominent signal at LHC and currently subject to intensive 
searches at Tevatron. The gluon-fusion channel (via a top-quark loop) is by far the largest produc- 
tion mode and known including the NNLO QCD corrections H54U65l - [67Tl . 

In Table [7] the total production cross sections for the Higgs-boson are presented as a function 
of the Higgs-boson mass m# at Tevatron and for a series of foreseen collision energies at the LHC 
(using m t = 173 GeV). The relative error from the PDF fit amounts to 5.5 - 10% at Tevatron and 
to 2.5 - 3% at the LHC at the higher energies and to 3.5 - 4.5% at yfs = 7 TeV. Again we do not 
consider the theoretical uncertainty due to the variation of the factorization and renormalization 
scale (typically amounting to ~ 9 - 10% at NNLO). In Figure [14] we compare the production cross 
sections to the results obtained using the PDFs of Ref. RTTl . The MSTW08 predictions yield higher 
values. For the LHC energies both analyses agree at lower Higgs masses Mh ~ 100 GeV and a 
gradual deviation reaching 3cr P at M# = 300 GeV of the MSTW08 values is observed. Our values 
at Tevatron are lower than those of MSTW08 by ~ 3o> in the whole mass range. At the LHC 
energies the difference can be attributed to different gluon PDFs and values for a s . The cross 
sections take very similar values for light Higgs masses, but beyond scales /j 2 ~ 10 4 GeV 2 the 
values obtained with MSTW08 are larger. 

6 Conclusions 

The precision of the DIS world data has reached a level which requires NNLO analyses to de- 
termine the PDFs and to measure the strong coupling constant a s (Mz). This also applies to the 
most prominent scattering processes at hadron colliders such as the Drell-Yan process, W ± -, Z- 
boson, Higgs-boson and top-quark pair-production. In the present analysis we have performed an 
NNLO fit to the DIS world data, Drell-Yan- and di-muon data along with a careful study of the 
heavy flavor effects in the DIS structure function F%. In the analysis we have taken into account 
correlated errors whenever available. In total, 25 parameters have been fitted yielding a positive 
semi-definite covariance matrix. With this information one may predict the error with respect to 
the PDFs, a s (Mz), mi, and m c for hard cross sections measured, including all correlations. For 
applications to hadron collider processes we have determined 3-, 4- and 5-flavor PDFs within the 
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ma 1 Ge V 


Tevatron 


LHC 7 TeV 


LHC 10 TeV 


LHC 14 TeV 


100. 


1.381 ±0.075 


21.19±0.58 


39.17+1.05 


67.28 ±1.77 


110. 


1.022 ±0.061 


17.30 ±0.49 


32.52 ±0.88 


56.59 + 1.51 


120. 


0.770 ±0.049 


14.34 ±0.41 


27.38 ±0.72 


48.25 ±1.23 


130. 


0.589 ±0.041 


12.03 ±0.36 


23.33 ±0.61 


41.60 ±1.07 


140. 


0.456 ±0.033 


10.21 ±0.31 


20.08 ±0.55 


36.23 ±0.92 


150. 


0.358 ±0.028 


8.75 ±0.27 


17.45 ±0.48 


31.83±0.82 


160. 


0.283 ± 0.024 


7.56 ±0.24 


15.29 ±0.43 


28.20 ±0.72 


170. 


0.226 ± 0.020 


6.59 ±0.21 


13.51 ±0.37 


25.16±0.65 


180. 


0.183 ±0.017 


5.78 ±0.19 


12.01 ±0.35 


22.60 ±0.60 


190. 


0.148 ±0.014 


5.11 ±0.17 


10.75 ±0.31 


20.44 ±0.53 


200. 


0.121 ±0.013 


4.55 ±0.16 


9.69 ±0.28 


18.59 ±0.49 


210. 




4.07 ±0.15 


8.78 ±0.26 


17.01 ±0.44 


220. 




3.67 ±0.14 


8.00 ±0.24 


15.64 ±0.42 


230. 




3.32±0.13 


7.33 ±0.22 


14.46 ±0.38 


240. 




3.02 ±0.12 


6.75 ±0.21 


13.44±0.37 


250. 




2.77 ±0.11 


6.25 ±0.20 


12.55 ±0.35 


260. 




2.55 ±0.10 


5.82±0.19 


11. 79 ±0.32 


270. 




2.36 ±0.10 


5.45 ±0.18 


11.12 + 0.31 


280. 




2.19±0.10 


5.13 + 0.17 


10.56 ±0.30 


290. 




2.06 ±0.09 


4.86 ±0.17 


10.08 ±0.29 


300. 




1.94 ±0.09 


4.63 ±0.16 


9.69 ±0.28 



Table 7: The total cross sections for Higgs-boson production [pb] at Tevatron and the LHC at the scale 
/u = Mh with the uncertainties estimated from the fit results in the present analysis. 



GMVFN scheme applying the BMSN description. We have performed a detailed study of the 
heavy flavor contributions to deep inelastic scattering comparing to experimental data. We have 
compared to different treatments used in the literature and found that both the FFN scheme and 
the BMSN scheme yield a concise description of the DIS data at least for the kinematic range 
of HERA, and that no modifications of these renormalization group-invariant prescriptions are 
needed. In the present analysis we have obtained a ? (M|) with an accuracy of ~ 1.5%. The values 
quoted in Eqs. d29l) and (|3Q|) are found to be in very good agreement with the non-singlet analysis 
of Ref. 114211 . which relied on a sub-set of the present data only, and with the results of Ref. [|45ll . 
The central value of a s (M 2 ) steadily converges going from LO to NLO to NNLO, or even to N 3 LO 
in the non-singlet case 11421 . The differences in the central values (determined at ju 2 = Q 2 ) provide 
a good estimate of the remaining theory errors. It is very hard to achieve a better accuracy on 
a s (M 2 ) than obtained at the moment, given the theoretical uncertainties (reaching values around 
~ 0.7%), which arise from the difference between the FFN and BMSN scheme, from quark mass 
effects, from 4-loop effects in the strong coupling constant from (the yet unknown) effect of the 
4-loop singlet anomalous dimensions, or from remainder higher twist effects and so on. However, 
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100 120 140 160 180 200 220 240 260 280 300 

M H (GeV) 

Figure 14: The lap error band for the Higgs-boson production cross sections [pb] at Tevatron and the 
LHC at the scale of = Mh employing the PDFs from the present analysis (shaded area) in comparison with 
the central values for the case of PDFs of Ref. H71 (dash-dotted lines). 

potential high-luminosity measurements planned at future facilities like EIC [68], requiring an ex- 
cellent control in the systematics, may provide future challenges to the precision on the theoretical 
side. 

We have discussed the NNLO PDFs of the present fit and compared to other global analyses. A 
comparison to the results of MSTW08 in the region fj 2 = 4 to 10 4 GeV 2 show that smaller values 
for the light PDFs for lower values of x are obtained in Ref. Il47l . Moreover, the gluon distribution 
of Ref. 11471 at low scales p? = 4 GeV 2 does strongly deviate from ours turning to negative values 
at x ~ 5 • 1(T 5 . At large values of x the gluon distribution of Ref. WT\ is slightly larger than ours. 
Somewhat smaller values are also obtained for the c- and £-quark distributions. The PDFs obtained 
in Ref. 11431 agree very well with the results of the present analysis. 

We have illustrated the implications of the PDFs for standard candle processes, such as W ± - 
and Z-boson production at hadron colliders. Comparison to MSTW08 yields a 2o> lower result 
for Tevatron and better agreement is obtained for the LHC energies. Conversely, the inclusive tt 
production cross section of both analyses agree at Tevatron energies, but for the LHC larger results 
by ~ Icrp are obtained with MSTW08. For the inclusive Higgs-boson production cross section at 
Tevatron the PDFs of MSTW08 yield a 3crp larger value in the whole mass range, while for LHC 
energies both predictions agree for masses M# ~ 100 GeV, and MSTW08 gives by 3o> larger 
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values for M# ~ 300 GeV. Of course, all observed differences have to be considered in view of the 
statistical and systematic accuracies finally to be obtained in the experimental measurements. 

The PDFs of the present analysis allow for detailed simulations of the different inclusive pro- 
cesses at the LHC and are of central importance in monitoring the luminosity. Precision mea- 
surements of inclusive processes at hadron colliders open up the opportunity to further refine the 
understanding of the PDFs of nucleons. This applies to both, the final analyses at Tevatron and the 
future measurements at the LHC. During the last years our understanding of PDFs has steadily im- 
proved at the NNLO level and upcoming high luminosity data from hadron colliders will continue 
in this direction. 

Grids, which allow fast access to our 3-, 4-, and 5-flavor PDFs in a wide range of x and Q 2 
(including the PDF uncertainties considered) are available online at |[69l . 
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